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Abstract 

The relationship between statics and dynamics proposed by Franz, 
Mezard, Parisi and Peliti (FMPP) for slowly relaxing systems [Phys. Rev. Lett. 
81, 1758 (1998)] is investigated in the framework of non disordered coarsen- 
ing systems. Separating the bulk from interface response we find that for 
statics to be retrievable from dynamics the interface contribution must be 
asymptotically negligible. How fast this happens depends on dimensional- 
ity. There exists a critical dimensionality above which the interface response 
vanishes like the interface density and below which it vanishes more slowly. 
At d = 1 the interface response does not vanish leading to the violation of 
the FMPP scheme. This behavior is explained in terms of the competition 

between curvature driven and field driven interface motion. 
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Recently Franz, Mezard, Parisi and Peliti (FMPP) have proposed a connection 
between static and dynamic properties of slowly relaxing systems. This is of fundamental 
importance for disordered systems such as spin glasses, since the low temperature equilibrium 
properties are hard to compute and still object of controversy after many years of intensive 
study. The existence of a bridge between statics and dynamics then offers a useful alternative 
tool for the investigation of the equilibrium state. 

Here we are interested in the analysis of this connection for non disordered coarsening 
systems, such as a ferromagnet quenched below the critical point. In this case the structure 
of the equilibrium state is simple and well known. Therefore, these systems are particularly 
suitable for a detailed understanding of the method. Nonetheless, existing results are quite 
puzzling. For the nearest neighbors Ising model numerical results in space dimension d > 2 
fit into the FMPP scheme, while recent exact analytical results in the d = 1 case ]3J show 
a qualitatively different behavior excluding any connection between the relaxation properties 
and the structure of the equilibrium state. In this paper we investigate the problem through 
a careful analysis of the linear response function as the space dimensionality is changed. This 
allows to put together the rich and interesting picture of what goes on in the off equilibrium 
response of a coarsening system and to uncover the mechanism whereby the FMPP scheme 
is or is not verified. 

Let us first outline the problem. We consider a phase ordering system with a scalar order 
parameter quenched at the time t = from high temperature to a final temperature T below 
the critical point. The time evolution takes place by formation and coarsening of domains 
of the opposite equilibrium phases 0. The characteristic feature of the process is the 
coexistence of fast and slow dynamics. Within domains local equilibrium is reached rapidly 
while the coarsening process is slow. In the case of non conserved order parameter, as it 
will be considered here, the typical domain size grows like L(t) ~ t 1 ^ 2 . This phenomenology 
suggests the split of the order parameter into the sum of two independent components |6| 




(1) 
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where ip(x,t) describes equilibrium thermal fluctuations within domains, the ordering com- 
ponent a(x, t) takes values im^ in the bulk of domains with the change of sign occurring at 
the interfaces and ttit is the equilibrium spontaneous magnetization. The split (|1|) accounts 
well for the observed behavior [0] of the autocorrelation function C(t,t') = (cf)(x,t)(f)(x,t')) 
where t > t' > 0. This quantity can be written as the sum of two terms 

C(t,t') = C st (t-t') + C ag (^j (2) 

representing, respectively, the stationary dynamics of thermal fluctuations ij)(x,t) and the 
slow out of equilibrium dynamics of a(x, t) obeying an aging form. At equal times C st (t—t' = 
0) = itlq — m^, where mo is the T = spontaneous magnetization, and C ag (l) = m^. 
Furthermore, due to the wide separation of time scales, in the range of time over which C st 
decays to zero the aging contribution remains practically constant C ag (t/t') ~ m^. 

Suppose, next, that at the time t w > a small random field with expectations h(x) = 0, 
h(x)h(y) = hl5{x — y) is switched on. Given the structure (|1|), the perturbation affects only 
the ordering component a(x, t), leaving thermal fluctuations unaltered, exactly as it happens 
in the equilibrium ordered phase under the action of an external field. Therefore, each 
unperturbed configuration is mapped into a new one (p(x,t) — > (ph(x,t) = ip(x,t) + ah(x,t) 
with a modified ordering component. In general, the perturbation will modify both the bulk 
and the interface behavior of a(x,t). Accordingly, to linear order we may write 

a h (x,t) = a(x,t) + J <Ix'xb{x — , t, t w )h{x') 

+ J ' dx'x^x - X* \t,t w )h(x') (3) 

where the bulk response function xb, accounting for the change in the magnetization within 
domains, must be related to thermal fluctuations via the equilibrium fluctuation dissipation 
theorem (FDT) 

Txb(x - x',t,t w ) = 

C st (x -x',t-t w = 0)- C st (x -x',t- t w ). (4) 
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At this point one can already get a glimpse on what conditions must be realized for statics 
and dynamics to be connected. Since the bulk response function xb probes the equilibrium 
fluctuations and that is where the information on the equilibrium state is stored, for the 
FMPP scheme to work the interface contribution \i must disappear. Current belief is 
that indeed this is what happens assuming that xi g° es like the interface density 

pi(t) ~ L~ l [t) ~ t~ x l 2 . However, to a closer scrutiny things are not so straightforward 
and the interface contribution turns out to have more structure than hitherto believed. In 
particular, there is an unexpected dependence on dimensionality which, in the end, accounts 
for the discrepancies in the Ising model mentioned above. Before entering the calculation of 
Xi, it is useful to complete the outline of the problem illustrating a few points. 

i) FMPP scheme. Averaging over the external field and the noise the on site total 
response function is given by 

1 
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X(t, t w ) = I dx((p h (x, t))h(x) = XB(t, t w ) + xi(t, t w ). (5) 



In the FMPP scheme one assumes that this quantity obeys the out of equilibrium general- 
ization of the FDT proposed by Cugliandolo and Kurchan |]], namely that for large times 
x(t,t w ) depends on the time arguments through the autocorrelation function x(C(t,t w )). 
Furthermore, one assumes that lim^oo t w ) = Xeg, where Xeg is the equilibrium response 
function. With these two hypothesis, the connection between statics and dynamics is derived 
in the form 

D(q) = P(q) (6) 

where D(q) = —T F d ^ c is the dynamical quantity, while P(q) = lim^o Ph{q) is the 
static one. Ph{q) is the overlap distribution in the perturbed Gibbs state. We recall that for 
a ferromagnetic system in the unperturbed Gibbs state below Tq the overlap distribution is 
given by |10 



P{q)= l -5{q-ml)+ l -5{q + ml). (7) 
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ii) Validity of the FMPP scheme. If in (||) we take into account only the bulk contribution 
Xb-i neglecting Xi, then the FMPP scheme is verified. In fact, using (0) and recalling that 
C ag (t/t') ~ m\ in the range of times where C s t(t — t') decays to zero, from (^) follows 

I C(t, t) - C(t, t w ) for m\<C <ml 
T X (C) = t o (g) 

TOq — m\ for C < m\. 

Eq. ([]) then implies P{q) = 5{q — m\). The presence of one ^-function instead of two as 
in (|?p does not mean that in the limit h — > symmetry is broken. Rather, it means that 
from a response due only to thermal fluctuations it is impossible to distinguish one pure 
state from the other and everything goes as if symmetry was broken. Numerical simulations 
of the Ising model for d > 2 show evidence for convergence toward the structure (|§D in 
the parametric plot of x versus C as t w becomes large (Fig. 1). This indicates that the 
interface term xi i n © must be asymptotically negligible. 

iii) Violation of the FMPP scheme. From the exact analytical solution of the d = 1 Ising 
model one finds a different behavior. In d — 1 it is important to realize that in order to make 
compatible the existence of a linear response regime, which requires ^ << 1 and therefore 
T > 0, with the existence of an equilibrium state of the form (0), one must take the limit 
J — > oo in the ferromagnetic coupling. This in turn implies m\ = m\. Hence, if (|8|) were to 
hold, Tx(t,t w ) would vanish. Instead, one finds j|] Tx(C) = (v^/tO arctan v^cot (ttC/2) 
which yields D(q) = it cos (nq/2) sin (nq/2) / 2 — sin (\q) and which is in no way related 



to (0), as shown in Fig. 1. Furthermore, notice that J = oo leads to the suppression 
of thermal fluctuations within domains. Therefore, from the above analytical form and 
x{t,t w ) = xi{t,t w ) follows lim^oo T X i(t,t w ) = lim c ._ >0 T X {C) = l/y/2. Therefore, the 
violation of the FMPP scheme in d = 1 is due to an asymptotically dominant interface 
contribution. 

Having made clear the necessity of investigating the relative importance of the bulk and 
interface terms as dimensionality is varied, we now introduce a semiphenomenological model 
for t w ) . This is based on the standard methods of the late stage theory in phase ordering 
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kinetics [[|, when only interface motion is of interest. Dropping the bulk term in (|^) and 
defining 07(2, t) = a(x, t) + J dx'xi(x — x 1 ', t, t w )h(x') we resort for this quantity to the time 
dependent Ginzburg-Landau model without thermal noise 

do" 

-^ = V 2 a I + m 2 T a I -a] + h(x). (9) 

Next, in order to allow for the action of the field on the interface motion, while keeping the 
domain saturation fixed at the unperturbed level ±mf, we make the ansatz 

° A ™ = -$&» (10) 



and we make an approximation of the gaussian auxiliary field (GAF) type |5,11|. The idea is 
that the non-linearity of the transformation (|T0|) is enough to take care of the non-linearity in 
the problem and the auxiliary field u(x, t) can be treated in mean field theory. Inserting ( |TT)| ) 
in and linearizing the equation of motion for u(x,t) we find 

^ = V 2 « + m^-^((Vtt) 2 )+/i(x) (11) 
at ' mi, 



where ((Vm) 2 ) is evaluated self-consistently. The key point is that the linear equation ( |TTD 
allows u(x, t) to grow unboundedly yielding via (|l~0|) 07(2, t) ~ rriTsign[u(x, t)} which enforces 
the saturation of domains at the required unperturbed value. Making a further mean field 
approximation through the replacement of ([T0|) by 07(2, t) = m T ^y= with S(t) = (u 2 (x, t)), 



we have xi(t,t w ) = rriTXu(t } t w ) / \J S(t) where Xu is the response of the auxiliary field. 
Computing Xu(t, t w ) and S(t) from ([11]) and defining Xeff{t, t w ) by Xi(t, t w ) = pi(t)x e ff{t, t w ) 
we find our main result 

Xeff(t,t w ) =t l ~iF{t w /t) (12) 



with F(x) = Afi dyy~—(l -y + ^)- d / 2 . Here A is a dimensionality dependent constant 
and to is a microscopic time related to the momentum cutoff A by to = A -2 . The meaning 
of Eq. (|l^) can be made transparent regarding the response of the system as due to a set of 
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interfaces each contributing through an effective response Xeff- From (|T2|) follows that this 
quantity obeys 

X eff(t,t w )~(t-t w ) a (13) 

with a = for d > 2, a = 1 - d/2 for d < 2 and Xeff ~ log(t - t w ) for d = 2 0. Eq. (|lj) 
applies both in the short (t — t w « t w ) and in the large (t — t w » t w ) time region, with 
a change in the prefactor taking place about t — t w ~ The full time dependence of 
Xeff(t,t w ) obtained by plotting out Eq. (|l^) for different values of d is displayed in Fig. 2a. 
A completely analogous behavior is obtained in the Ising model. We have computed Xi(t, t w ) 
for d = 1,2,3,4 in the Ising case by suppressing spin flips in the bulk of domains and we 
have plotted Xeff(t, t w ) in Fig. 2b [0 . The common features of Fig. 2a and Fig. 2b may be 
summarized stating that in both cases Xeff obeys the power law ([U|) and that there exists a 
critical dimensionality d c such that the exponent a is zero for d > d c while it grows positive 
with decreasing dimensionality for d < d c , reaching in both cases the final value a = 1/2 
at d = 1. The differences are that d c = 2 in the GAF approximation, while from the Ising 
simulations there is a good evidence (Fig. 2b) for d c = 3. One particular consequence of 
this, on which we comment below, is that while in the former case we have a = with 
logarithmic growth at d = 2, instead in the d = 2 Ising model one finds a = 1/4 (Fig. 2b). 

The interpretation of the above results goes as follows. In the perturbed system interface 
motion is driven not only by curvature, as in the unperturbed case [§], but also by the 
external field. These two mechanisms compete. For d > d c the curvature mechanism 
dominates. The external field then affects only the spins strictly belonging to the interface. 
Within a microscopic time Xeff saturates (see lines for d = 3 in Fig. 2a and for d = 4 in 
Fig. 2b) leaving the overall response xi t° decrease as pi(t). Instead, if dimensionality is 
low enough (d < d c ) to weaken the curvature mechanism to the point that the field driven 
motion may start to play a role in the response, then the single interface response Xeff 
grows with time like t a counteracting the default decrease in xi due to pi(t). This effect 
can be understood realizing that when interfaces are field driven the field produces a large 



scale optimization of domain positions with respect to field configurations fll4| . d = 1 is 
the extreme case where there is no more curvature mechanism. Then, interface motion is 
entirely field driven yielding Xeff which grows like t 1 ^ 2 and compensates exactly for piit) 
producing a non vanishing lim^oo Xi(ty tw)- 

In summary, we have investigated the relationship between the out of equilibrium re- 
sponse function and the structure of the equilibrium state for coarsening systems. For the 
FMPP scheme to hold, the out of equilibrium response must be dominated by equilibrium 
fluctuations in the bulk of domains. This we have made precise by separating the bulk from 
interface response and by analyzing the behavior of the interface contribution xi{t,tw)- On 
the basis of a GAF model we have found that there exists a critical dimensionality d c = 2 
such that for d > d c Xii^i^w) behaves like the interface density pi(t), while for d < d c it 
vanishes slower than pi(t) and does not vanish anymore at d = 1, yielding the violation of 
the FMPP scheme. We have explained this behavior identifying d c as the dimensionality 
below which the external field competes effectively with the curvature in driving interface 
motion. The overall picture is confirmed by numerical results for the Ising model, apart 
from the upward shift from d c = 2 to d c = 3 in the critical dimensionality. This means that 
in the Ising case the field driven mechanism competes with the curvature even more effi- 
ciently than in the GAF approximation. In particular, the field driven mechanism is clearly 
observable in the d = 2 Ising model where a = 1/4 makes the interface response much 
less preasymptotic than what was estimated on the basis of the interface density argument. 
Finally, the counterpart of the statement that a persistent interface contribution leads to 
the violation of the FMPP scheme, is that even if Xi{t>tw) vanishes asymptotically, due to 
the field driven interface motion it may vanish so slowly to hide the realization of the FMPP 
scheme. 
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FIG. 1. Tx(t,t w ) for the Ising model. For d = 2 the curve shows the expected behavior on the 
basis of Eq. (|8|) for t w —* oo with T = 2.2, J = 1. For ci = 1 the curve is the plot of the analytic 
form of Ref. [4]. The inset shows the corresponding D(q). 
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a) GAF model 





FIG. 2. Xeff{t, tw) hi the semiphenomenological model with t w = 10 8 (a) and in the Ising model 
without spin flips in the bulk (b). For d = 1, 2, 3, 4, the temperature, waiting time and linear system 
sizes L of the simulations are T = 0.48, 2.2, 3.3, 4.4, t w = 10 3 , 10 3 , 10 2 , 10 and L = 10 6 , 512, 128, 42 
with J = 1 and averages over 170, 6045, 2590, 922 realizations. The dashed lines are power laws with 
the corresponding exponent a. For d = 3 in (b) the curve is very well fitted by 0.33+0.066 \og{t—t w ). 
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